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ABSTRACT 

Granato et al. (2004a) have elaborated a physically grounded model exploiting the 
mutual feedback between star-forming spheroidal galaxies and the active nuclei grow- 
ing in their cores to overcome, in the framework of the hierarchical clustering scenario 
for galaxy formation, one of the main challenges facing such scenario, i.e. the fact 
that massive spheroidal galaxies appear to have formed much earlier and faster than 
predicted by previous semi-analytical models, while the formation process was slower 
for less massive objects. After having assessed the values of the two parameters that 
control the effect of the complex and poorly understood radiative transfer processes 
on the time-dependent spectral energy distributions (SEDs), we have compared the 
model predictions with a variety of IR to mm data. Our results support a rather strict 
continuity between objects where stars formed (detected by (sub)-mm surveys) and 
evolved massive early-type galaxies, indicating that large spheroidal galaxies formed 
most of their stars when they were already assembled as single objects. The model 
is remarkably successful in reproducing the observed redshift distribution of ^ 20 
galaxies at 2; > 1, in contrast with both the classical monolithic models (which over- 
estimate the density at high-z) and the semi-analytic models (that are systematically 
low) , as well as the ratio of star- forming to passively evolving spheroids and the counts 
and redshift distributions of Extremely Red Objects (EROs), although the need of a 
more detailed modelling of the star formation history and of dust geometry is indi- 
cated by the data. The model also favourably compares with the ISOCAM 6.7 /im 
counts, with the corresponding redshift distribution, and with Spitzer/IRAC counts, 
which probe primarily the passive evolution phase, and with the (sub)-mm SCUBA 
and MAMBO data, probing the active star-formation phase. The observed fraction 
of 24 /im selected sources with no detectable emission in either the 8 /im or i?-band 
(Yan et al. 2004a) nicely corresponds to the predicted surface density of star-forming 
spheroids with 8 /im fluxes below the detection limit. Finally, distinctive predictions 
for the redshift distributions of 24 /im sources detected by Spitzer/MIPS surveys are 
pointed out. 

Key v^rords: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
formation - QSOs: formation 



1 INTRODUCTION 

The standard Lambda Cold Dark Matter (ACDM) cosmol- 
ogy is a well established framework to understand the hi- 
erarchical assembly of dark matter (DM) halos. Indeed, it 
has been remarkably successful in matching the observed 
large-scale structure. However the complex evolution of the 
baryonic matter within the potential wells determined by 
DM halos is still an open issue, both on theoretical and on 
observational grounds. 

Full simulations of galaxy formation in a cosmological 



setting are far beyond present day computational possibili- 
ties. Thus, it is necessary to introduce at some level rough 
parametric prescriptions to deal with the physics of baryons, 
based on sometimes debatable assumptions (e.g. Binney 
2004). A class of such models, known as semi-analytic mod- 
els, has been extensively compared with the available in- 
formation on galaxy populations at various redshifts (e.g. 
Lacey et al. 1993; Kauffmann, White & Guiderdoni, 1993; 
Cole et al. 1994; Kauffmann et al. 1999; Somerville & Pri- 
mack 1999; Cole et al. 2000; Granato et al. 2000; Benson et 
al. 2003; Baugh et al. 2004). 
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Figure 1. Contributions of spheroidal galaxies to the K-haad number counts and redshift distributions of galaxies brighter than K = 20, 
23, 24. The dashed, solid and dotted lines refer to fc = 3, 30, 300 M^''^ respectively (see Sect.|^. We adopt k = 30Mg^^ as our reference 
value. Data for K band counts are from Moustakas et al. (1997), Kochanek et al. (2001), Saracco et al. (2001), Totani et al. (2001), 
Cimatti et al. (2002). In the upper-right panel (the K< 20 redshift distribution), the thin solid and the thin solid with crosses histograms 
(the redshift bin Az is 0.15) are respectively from Cimatti et al. (2002) and Somerville et al. (2004) (scaled to the same area of the K20 
survey). In the lower panels the thin continuous histograms (A2=0.3) arc from Kashikawa et al. (2003). 



The general strategy consists in using a subset of ob- 
servations to calibrate the many model parameters provid- 
ing a heuristic description of baryonic processes we don't 
properly understand. Besides encouraging successes, current 
semi-analytic models have met critical inconsistencies which 
seems to be deeply linked to the standard recipes and as- 
sumptions. These problems are in general related to the 
properties of elliptical galaxies, such as the color-magnitude 
and the [a/Fe]-M relations (Cole et al. 2000; Thomas 1999; 
Thomas et al. 2002), and the statistics of sub- mm and deep 
IR selected (I- and K-band) samples (Silva 1999; Chapman 
et al. 2003; Kaviani et al. 2003; Daddi et al. 2004; Kashikawa 
et al. 2003; Poh et al. 2003; Pozzetti et al. 2003; Somerville 
et al. 2004). Discrepancies between the observationally es- 
timated merger rate evolution and that predicted by semi- 
analytic models have also been reported (Conselice et al. 
2003). 

However, the general agreement of a broad variety of ob- 
servational data with the hierarchical scenario and the fact 



that the observed number of luminous high-redshift galaxies, 
while substantially higher than predicted by standard semi- 
analytic models, is nevertheless consistent with the number 
of sufficiently massive dark matter halos, indicates that we 
may not need alternative scenarios, but just some new in- 
gredients or assumptions (see e.g. Baugh et al. 2004; Tecza 
et al. 2004). 

Previous work by our group (Granato et al. 2001; Ro- 
mano et al. 2002; Granato et al. 2004a) suggests that a cru- 
cial ingredient is the mutual feedback between spheroidal 
galaxies and active nuclei at their centers. Granato et al. 
(2004a, henceforth GDS04) presented a detailed physically 
motivated model for the early co-evolution of the two com- 
ponents, in the framework of the ACDM cosmology. 

In this paper, we present a comparison of the model 
with the observed number counts and redshift distributions 
in several near-IR (NIR) to (sub)-mm bands. 

In Section |5| we give a short overview of the GDS04 
model for spheroidal galaxies, and we discuss the determina- 
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tion of the parameters controlling the time-dependent spec- 
tral energy distributions (SEDs) of spheroidal galaxies. In 
Section El > the model predictions are compared with a va- 
riety of recent data. The main conclusions are summarized 
and discussed in Section 2] 

We adopt the "concordance" cosmological model, with 

= 0.3, = 0.7, Ho = 70kms"^Mpc~\ 



2 MODEL DESCRIPTION AND ASSESSMENT 
OF THE PARAMETERS 

2.1 The GDS04 model 

While referring to GDS04 for a full account of the model as- 
sumptions and their physical justification, we provide here, 
for the reader's convenience, a brief summary of its main 
features. 

The model follows with simple, physically grounded 
recipes and a semi-analytic technique the evolution of the 
baryonic component of proto-spheroidal galaxies within 
massive dark matter (DM) halos forming at the rate 
predicted by the standard hierarchical clustering scenario 
within a ACDM cosmology. The main novelty with respect 
to other semi-analytic models is the central role attributed to 
the mutual feedback between star formation and growth of a 
super massive black hole (SMBH) in the galaxy center. Fur- 
ther relevant differences are the assumption that the large 
scale angular momentum does not effectively slow down the 
collapse of the gas and the star formation in massive halos 
virialized at high redshift, and the allowance for a clump- 
ing factor, C ~ 20, speeding up the radiative cooling so 
that even in very massive halos (Mvir ~ 10^^ Mq), the gas, 
heated to the virial temperature, can cool on a relatively 
short (~ 0.5-1 Gyr) timescale, at least in the dense central 
regions. 

The GDS04 prescriptions are assumed to apply to DM 
halos virializing at Zvii k, 1.5 and Mvir ^ 4 x lO^^A/0. These 
cuts are meant to crudely single out galactic halos associ- 
ated with spheroidal galaxies. Disk (and irregular) galaxies 
are envisaged as associated primarily to halos virializing at 
Zvii £ 1.5, some of which have incorporated most halos less 
massive than 4 x 10^^ Mq virializing at earlier times, that 
may become the bulges of late type galaxies. However, the 
model does not address the formation of disk (and irreg- 
ular) galaxies and these objects will not be considered in 
this paper. On the other hand, we do not expect that proto- 
spheroidal galaxies in the star-forming phase have already a 
spheroidal shape; rather they are likely to have an irregular, 
complex, disturbed morphology. 

The kinetic energy fed by supernovae is increasingly ef- 
fective, with decreasing halo mass, in slowing down (and 
eventually halting) both the star formation and the gas ac- 
cretion onto the central black hole. On the contrary, star 
formation and black hole growth proceed very effectively in 
the more massive halos, giving rise to the bright SCUBA 
phase, until the energy injected by the active nucleus in the 
surrounding interstellar gas unbinds it, thus halting both the 
star formation and the black hole growth (and establishing 
the observed relationship between black hole mass and stel- 
lar velocity dispersion or halo mass). Not only the black 
hole growth is faster in more massive halos, but also the 
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Figure 2. Contribution of star forming (dot-dashed line) and 
passively evolving (long-dashed line) spheroids to the K < 20 
redshift distribution (the same as in Fig.0. The thick solid line 
is the model total, the thin line the Cimatti et al. (2002) data, 
and the thin line with crosses the data from Somerville et al. 
(2004) scaled to the same area. Note that Somerville et al. give 
magnitudes in the AB system (-R'as = K + 1.85). 



feedback of the active nucleus on the interstellar medium is 
stronger, to the effect of sweeping out such medium earlier, 
thus causing a shorter duration of the active star-formation 
phase. As a result, in keeping with the previous proposition 
by Granato et al. (2001), the physical processes acting on 
baryons reverse the order of formation of spheroidal galax- 
ies with respect to the hierarchical assembling of DM halos, 
leading to the Anti-hierarchical Baryon Collapse (ABC) sce- 
nario. 



2.2 Spectral energy distribution of spheroidal 
galaxies 

As in GDS04, we compute the SEDs of spheroidal galaxies 
using our code GRASIL, described in Silva et al. (1998)^. 
GRASIL computes the time-dependent UV to radio SED 
of galaxies, given their star formation and chemical enrich- 
ment histories (derived as described in GDS04), and with 
a state of the art treatment of dust reprocessing. The lat- 
ter point is fundamental, since during the phase of intense 
star-formation, young stars are mixed with a huge amount 
of gas, quickly chemically enriched and dust polluted. The 
use of local templates (such as M82 or Arp220) may be in- 
appropriate in these extreme conditions. 

One of the most important distinctive features of 
GRASIL is that it has included, for the first time, the effect 
of differential extinction of stellar populations (younger stel- 
lar generations are more affected by dust obscuration), due 



^ The code is av ailable at \http://adlibitum.oat.ts. astro.it/silvaj 
default.html or a.i\http://web.pd. astro. it/granato\ 
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Figure 3. Possible contributions of spheroids to the observed 15 and 850 lira number counts and redshift distributions (Sis^m > O.lmjy 
within 6 X 10"'^ deg^ and Az = 0.2; Sgso/jm > Imjy). The dotted, sohd and dashed hnes refer to T)y[c(lAim) — 10, 60, 120 respectively. 
The 15 /xm data are from Elbaz et al. (1999), Gruppioni et al. (2002), Elbaz et al. (2002, thin histogram with error bars), Franceschini 
et al. (2003, thin histogram), Metcalfe et al. (2003). The 850 /xm data are from Hughes et al. (1998), Barger, Cowie, & Sanders (1999), 
Blain et al. (1999), Eales et al. (2000), Borys et al. (2002), Chapman et al. (2002). A low tmc(1A'™) implies an excessive contribution 
from spheroids to the 15 /xm counts, which are dominated by late-type, mainly starburst, galaxies at z ^ 1 (Hammer et al. 2004). On 

the other hand, a very high tmc(1/^"^) overpredicts the 850^m counts. Our reference model has tmc (l/^™) = 60. See text for details. 



to the fact that stars form in a denser than average environ- 
ment, the molecular clouds (MCs), and progressively get rid 
of them. 

Radiative transfer results depend on the poorly known 
geometry of the system. Broadly speaking, predictions of 
fluxes during the active star forming phase become more and 
more affected by the ensuing model uncertainties at shorter 
and shorter wavelengths. In GDS04 we presented predictions 
for the dust-enshrouded SCUBA galaxies, or for the essen- 
tially dust-free passively evolving galaxies; in both cases, the 
results were insensitive to the details of the dust distribution. 
On the other hand, key information on the early evolution 
of spheroidal galaxies is provided by near-IR (particularly 
Tf-band) and mid-IR (ISO and Spitzer) data. To test the 
model against them we need a more careful investigation of 
the dependence of our results on the structure of galaxies. 

In general, the GRASIL SEDs for spheroidal systems 



depend on five structural parameters: the fraction /mc of 
gas in the form of MCs rather than in the diffuse ISM (cir- 
rus), the optical depth of MCs for a source at their center 
T"MC (conventionally given at 1 /im hereafter), the escape 
timescale of newly born stars from MCs, t^, and the core 
radii of the distributions of diffuse dust, r^, and of stars 
outside MCs, r^, assumed to follow a King (1972) law. In 
the following we will take rj: = r*, and we will drop the 
superscript. 

However, in the conditions envisaged here and for the 
purposes of the present work, the results are significantly 
affected only by tmc- This is proportional to the dust to gas 
mass ratio 5 and to Mmc/Rmcj ^mc and -Rmc being the 
typical mass and radius of molecular clouds, respectively. As 
for 5, we have simply assumed that it increases linearly with 
the gas metallicity (e.g., Dwek 1998) and is ~ 1/100 (the 
Galactic value) for solar metallicity. Since the chemical evo- 
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Figure 4. Integral counts of star-forming spheroidal galaxies 
at 1200/im predicted by our reference model, compared with 
MAMBO data from Greve et al. (2004). 



lution of the galaxy is computed self-consistently, the true 
parameter is Mmc / ^mc > prefer to make reference to 

the value of tmc for solar metallicity (and thus 5 ~ 1/100). 
We have checked that all the results presented here, except 
those on Extremely Red Objects (EROs, see Sect. 13. 2t . are 
very little affected by variations, within rather wide ranges, 
of the other GRASIL parameters. 

However, an additional parameter is required to circum- 
vent the problem that, for stars already outside their parent 
molecular clouds, GRASIL takes into account only the rela- 
tively small optical depth due to cirrus, neglecting that due 
to other MCs present along the line of sight. This is not nec- 
essarily accurate in the extreme conditions occurring during 
the fast star forming phase of spheroids, particularly in the 
dense central regions. In general, the mean optical depth to a 
star outside its parent MC is Toutsidc = Tdrrus + (4/3) TV tmc, 
where we have assumed a uniform dust distribution inside 
MCs, and N is mean number of MCs along its line of sight. 



N = Sf 



M'gas,MC -RmC 

Mmc ri 



oc M^^sMcM^.^^^ , (1) 



1/3 

where Sf is a shape factor and we have assumed oc M^/. , 
Mvir being the virial mass of the galaxy. Whenever iV << 1, 
the optical depth drops from tmc to Tdrrus 1 at IR wave- 
lengths) when the stars move out of their parent MCs. This 
condition is usually met in nearby galaxies. On the contrary, 
in forming spheroids Toutsidc can be substantially larger than 
Tcirrus. The shapc factor Sf depends on the relative distribu- 
tion of stars and MCs, and its value is also a function of the 
age of the stellar population under consideration: younger 
stars (but old enough to be out of molecular clouds) are on 
average closer to regions rich in MCs. 

Thus, an accurate calculation of N is extremely diffi- 
cult, and we are forced to parameterize our ignorance. Also, 
the present release of GRASIL does not allow us to deal with 
the radiation transfer through the distribution of MCs. On 
the other hand, such detailed treatment is not crucial for our 
present purposes since whenever Toutsidc ^ tmc essentially 
all the starlight is reprocessed by dust. We have crudely ap- 
proximated Toutsidc as a step function: Toutsidc = Tcirrus for 
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Figure 5. 6.7^m counts and redshift distribution of sources with 
56.7 > 10/.tJy within 16 arcmin^. The thick solid line shows the 

1/3 

prediction of our reference model (i.e. fc = 30 ). The contri- 

butions of active and passive spheroids are represented by the dot- 
dashed and long-dashed lines respectively. The thin histogram in 
the lower panel shows the redshift distribution (with Az = 0.2) 
observationally estimated by Sato et al. (2004). The observed 
counts are from Metcalfe et al. (2003, squares) and Sato et al. 
(2003, diamonds). 



Mgas.MC ^ kM^h S'lld ^"outside = Tcirrus "I" TMC Othcrwise, 

treating k as an adjustable parameter. In order to take into 
account the dust heating by the absorbed starlight, in prac- 
tice we have kept the stars within their MCs for the full dura- 

2/3 

tion of the starburst, for galaxies where Mgas,MC > kMJ^^ . 

The parameter k affects mostly the results at A 10/^m. 
The A'-band counts and redshift distributions are best re- 
produced setting k ~ 30Mg'''^ (Fig-0- Note that spheroids 
dominate the redshift distributions at substantial redshifts, 
while leaving room for other populations at low redshifts. 
The figure shows also the effect of varying k. In particular, 
a low value of fc (fc ~ 10, not shown), would produce an 
almost perfect match of the K <2Q redshift distribution for 
z ^ 1, but the contributions at 7^ < 23 and K < 24: would 
decrease. A more detailed treatment of the complex behav- 
ior of the optical depth to stars outside MCs could improve 
the overall match, and will be considered for future releases 
of GRASIL. 

The results in Fig. Q refer to te = 0.01 Gyr, /mc = 
0.5 and Tc = 0.1(Mvir/10^^MQ)^/^ kpc, but change only 
marginally if te and /mc vary within rather broad inter- 
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Figure 6. if-band counts and redshift distributions for EROs {R — K > b). The thick solid line is the total predicted by our model, 
while the contributions of active and passive objects are shown by dot-dashed and long dashed lines respectively. Data from Daddi et 
al. (2000), Miyazaki ct al. (2003), Roche et al. (2002, 2003), Cimatti et al. (2003). The redshift histograms are for = 0.1 for K < 20, 
and for Az = 0.2 in the two lower panels. 



vals (5Myr < te < 0.5 Gyr, 0.2 < /mc ^ 0.9), and is in- 
creased or decreased by a factor of a few. 

As for TMC, Silva et al. (1998) found that the SEDs 
of starburst galaxies are well reproduced by GRASIL with 
TMC ~ a few tens (a figure in reasonable agreement with di- 
rect observations of Giant Molecular Clouds in the Galaxy, 
for which Mmc ~ W^Mq, tmc ~ 15pc). However, the phys- 
ical conditions in high-z star forming spheroids may be very 
different from those encountered in local systems, includ- 
ing starbursts and Ultra-Luminous IR Galaxies (ULIRGs). 
Thus, we have explored the generous range 3 < tmc < 200. 
This parameter controls the contribution (if any) of star- 
forming spheroidal galaxies to the bump in the 15 ^m ISO- 
CAM counts below a few mjy. Although the observed counts 
are affected by a considerable uncertainty, as shown by the 
spread of results from different surveys, and the modelling 
is complicated by the presence of PAH bands, they have 
been interpreted by phenomenological models in terms of 
evolving starbursts, with a minor contribution from spirals, 
particularly at the brighter flux density levels (e.g. Frances- 
chini et al. 2001; Gruppioni et al. 2002; Lagache et al. 2003; 



Hammer et al. 2004). Thus, the contribution from star- 
forming spheroids must be substantially suppressed, imply- 
ing TMC ~ 40 at 1/im (see Fig. EJ. A low value of tmc, and, 
correspondingly, a large contribution from spheroids to the 
sub-mjy bump of the 15 fim counts (upper- left panel), would 
entail a high-z tail of the redshift distribution of sources with 
Si5 > 0.1 mJy substantially larger than indicated by current 
data (lower left-hand panel). 

On the other hand, as illustrated by the right-hand pan- 
els of Fig. 13 an exceedingly large value of tmc would over- 
produce the 850 fim counts, while a very small value would 
underproduce them. The ensuing lower limit to tmc is how- 
ever model dependent, since the effect of decreasing it may 
be compensated by a smaller value of the dust emissivity 
index, that we take equal to 2. 

In the following we set tmc — 60 at 1/xm, a value yield- 
ing 850 /xm counts and redshift distributions consistent with 
the data, and a small, but non-negligible, contribution to 
the sub-mJy 15 ^m counts. 
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Figure 7. Dependence of predicted redshift distribution of EROs 
in the K20 survey on details of the model (see Fig.|£] upper right 
panel). The dashed line with crosses is obtained adding to the 
reference model (dashed line without crosses) a mild SF activity 
after the QSO phase (i.e. at the beginning of the passive evolution 
phase), at a level of 2% of the peak value and lasting 1 Gyr. The 
dot-dashed lines (for clarity these have been shifted by 0.5 in 
redshift) refer to objects before the QSO phase (i.e. in the active 
star- forming phase), adopting te = 10 Myr (line without crosses), 
or te = 50 Myr (with crosses). 



3 COMPARISON WITH THE DATA 

3.1 Counts and redshift distributions 

As discussed by Granato et al. (2001, 2004a), in the present 
framework the SCUBA-selected galaxies are interpreted as 
being mostly star-forming massive proto-spheroidal galaxies 
(Mvir > lO" *^ Mq). The values of the GRASIL parame- 
ters derived above nicely fit the 850 /xm SCUBA counts (see 
Fig- El right-hand panels), and the latest MAMBO 1200^m 
data (Fig. Ell. The model 850 fim counts are almost undistin- 
guishable from those already reported in Fig. 12 of GDS04, 
corresponding to a somewhat different choice of the pa- 
rameters, and the redshift distributions are also essentially 
unchanged. Indeed, the predictions in the sub-mm region 
are not very sensitive to the precise values of GRASIL pa- 
rameters, as far as proto-spheroidal galaxies are optically 
thick in the optical-UV. The predicted redshift distribu- 
tion of sources brighter than 1 mjy at 850 /^m has a broad 
peak at z ~ 2 with an extended tail towards higher red- 
shifts, due to star-forming spheroidal galaxies. As already 
remarked in GDS04, this redshift distribution is in good 
agreement with the limited information available so far. In 
particular, the median redshift of the Chapman et al. (2003) 
sample (5*850^111 ^ 5 mJy) is 2.4, with a quartile range of 
2 = 1.9 — 2.8. The median and quartile range for our model 
are 2.3 and 1.8 - 3.0. 

While the sub-mm counts are informative on the star- 
formation phase of the evolution of spheroidal galaxies, the 
7^-band counts (Fig.0| in the range 14 ^ K ^ 17 are domi- 
nated by the passive evolution phase. Fainter than if ~ 17 
also star-forming spheroids begin to show-up and become in- 
creasingly important at fainter magnitudes. The brighter (in 



terms of apparent magnitude) such objects are those of mass 
lower than the typical mass probed by present day sub-mm 
observations, i.e. those having, in the GDS04 model, a pro- 
longed star-formation phase extending to z ~ 1, and whose 
stars have negligible obscuration as soon as they come out 
of their parent MC (see Sect. 12.21 . This is illustrated by 
Fig. |5| showing that, for K ^ 20, the contribution from 
star- forming spheroids peaks at z 1. In particular, the 
predicted number of spheroids with ongoing star forma- 
tion at 1.4 < z < 2 and if < 20 is about twice that of 
passive spheroids, and their typical stellar mass is between 
5 X IO^^Mq and 3 x 10^^ Mq , in nice agreement with the find- 
ings of Daddi et al. (2004; see their Figure 5) and Somerville 
et al. (2004). The Star Formation Rate (SFR) of these ac- 
tive objects is ;$ lOOMQyr"^: according to our model, the 

K band selects spheroids when their SF (and dust obscura- 
tion) is declining or after it has stopped. Daddi et al. (2004) 
measured a dust-corrected average SFR of ~ 2OOM0yr~^, 
higher than our expectations, but still well below the peak 
SFR of spheroids according to our scenario. This may sug- 
gest that the present version of our model envisages a too 
drastic decline of SF activity in intermediate mass objects 
(see also Sect. 13.21 . Note also that Daddi et al. adopt a 
Salpeter IMF (oc Ai'^'^'^), while in the GDS04 model we 
have used a different IMF (cx M"^'^^ above 1 Mq, cx M"° * 
below) suggested by chemical constraints (see also Romano 
et al. 2002). This IMF implies a higher L/M ratio and, cor- 
respondingly, a SFR lower by about 30%. 

We predict that passively evolving spheroids almost dis- 
appear at if ^ 23. The observed high-z tails of the distri- 
butions for if ^ 23 and for if ^ 24 are thus mostly due to 
star-forming spheroids, but the counts are likely dominated 
by late-type galaxies, which could account for the big peak 
of the estimated (from photometric redshifts) 2 distribution 
in the range 0.2 ^ z ^ 0.9, as well as for the low-z shoulder 
at if < 20. 

As shown by Fig. |S] according to our model, spheroids 
may fully account for the observed 6.7/im ISO CAM counts 
and for the redshift distribution derived by Sato et al. (2004; 
see also Flores et al. 1999) . Above ~ 100 piJy they are pre- 
dicted to be mostly in the passive evolution phase, but the 
fraction of active star-forming spheroids increases with de- 
creasing flux, reaching ~ 50% at ~ 10 nJy. 



3.2 Extremely Red Objects 

Extremely red objects (EROs), with R ~ K > 5, have re- 
ceived special attention in recent years (Daddi et al. 2000; 
Smith et al. 2002; Roche et al. 2002, 2003; Cimatti et al. 
2002, 2003; Takata et al. 2003; Yan & Thompson 2003; Yan 
et al. 2004b; Webb et al. 2004; Moustakas et al. 2004), since 
their properties set crucial constraints on the early evolu- 
tionary phases of massive spheroidal galaxies. EROs are 
actually a mix of dusty star-forming and evolved galaxies 
formed at high redshifts (Cimatti et al. 2002, 2003; Mohan 
et al. 2002). 

We have worked out the counts and the redshift 
distributions of both star-forming and passively evolving 
spheroidal galaxies brighter than if = 19.2, 20, 20.3 that, 
according to our model, have ERO colours (Fig. |SJ. This 
prediction is the only one, among those presented in this 
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Figure 8. Contributions to the 1 — 1000/im background from 
passively evolving plus star-forming spheroidal galaxies (solid 
line). Passively evolving spheroids (dashed line) are important 
at A S: 15 /xm. Data are from Hauser & Dwek (2001) and Metcalfe 
et al. (2003). 



paper, which depends strongly also on GRASIL parameters 
other than tmc and k (see Sect. I2.2|l . Indeed, the predic- 
tions at 2 ^ 1 for dusty and star forming EROs require the 

computation of rest frame fluxes at A ^ 3000 A, which are 
severely dependent on details of dust obscuration. 

In sharp contrast with the standard semi-analytical 
models, which severely under- predict this population (e.g., 
Somerville et al. 2004), our model somewhat over-produces 
the number of EROs (of both classes), in the range 1.5 S 
2 $ 2. On the observational side, the present redshift dis- 
tributions may be biased by the strong clustering of these 
sources and affected by sampling variance, and some may be 
incomplete at z > 1 (Yan et al. 2004b). From the theoretical 
point of view, it should be noted that a mild star formation 
activity (say at a level ~ a few % of the peak rate), either 
due to residual gas or induced occasionally by interactions, 
neglected by our simple model, could significantly decrease 
the predicted number of passive galaxies with ERO colours, 
particularly at high redshifts. To illustrate this possibility. 
Fig. [7| shows the effect on predicted z-distribution of EROs 
brighter than isT = 20 of an 'artificial' (i.e. not predicted by 
the model equations) episode of star formation starting af- 
ter the feedback from the QSO has swept out the ISM (QSO 
phase in the terminology of GDS04), and lasting 1 Gyr. A 
decline of the SFR somewhat less abrupt than implied by 
the GDS04 model seems also to be suggested by the results 
of Daddi et al. (2004), who find that star forming objects at 
1.4 < < 2 in the K20 survey have a dust-corrected star for- 
mation somewhat higher than our prediction (Sect. l^TTl . As 
for dusty and star forming EROs, we reminded above that 
our prediction depends on details of dust obscuration. For 
instance a shorter and/or larger values of the galaxy scale 
length Tc would translate in less reddened objects, without 
affecting significantly the other results presented in this pa- 
per. An example of this is also shown in Fig. |7| 
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Figure 9. Upper panel: contributions of spheroidal galaxies to 
the 24/im counts predicted by our reference model compared 
with the Spitzer/MIPS data (Papovich et al. 2004). The long- 
dashed line shows the counts of passively evolving spheroids. 
Lower panel: predicted redshift distributions for S'24/jm ^ 0.45 
mjy (the SWIRE limit, Lonsdale et al. 2003), 0.1 and 0.01 mjy 



In conclusion, it is encouraging that our model has 
not troubles in producing enough high-z objects with ERO 
colours. More statistic in the observations will allow us to 
improve the treatment of some aspects of the model, such 
as the amount of star formation after the QSO phase and 
the dust obscuration. 



3.3 The IR background 

A further constraint on the model comes from the IR back- 
ground. In Fig. |H|we compare the 1 — 1000//m background 
spectrum yielded by our model with the data collected by 
Hauser & Dwek (2001). Longward of A ~ 200 /jm star- 
forming spheroids alone may explain the observed back- 
ground, while at shorter wavelengths, and, in particular, at 
the intensity peak, late-type (especially starburst) galaxies 
are expected to dominate. 
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that only few of these sources have Ss^m > 20 ^Jy. We ex- 
pect that the same is true at 0.7 ^m, although predictions 
at this short wavelength, probing the rest frame UV for our 
model sources, are very uncertain. 

In Fig. 1101 we have also compared the redshift distribu- 
tion predicted by Chary e al. (2004, their Fig. 4) for 5*24^1X1 > 
0.02 mjy with the one from our model for spheroids. Chary 
et al. take into account low luminosity starbursts, LIRGs and 
ULIRGs. Note that our model leaves room for (but does not 
include) populations other than forming spheroids only at 
redshift ^ 1.5; therefore we predict much fewer sources at 
high redshift as compared to Chary et al., but with a tail 
extending to higher z. Also according to our scenario, the 
populations responsible for the bulk of the 15 and 24 /xm 
counts should be at redshift ^ 1.5, since all massive halos 
virialized at higher z are used to build spheroids (see GDS04 
for more details). The observed 24 /xm redshift distribution 
will be a strong test for our model. 



Figure 10. Redshift distribution for S24fim > 0.02 mJy. The 
thick histogram corresponds to our reference model for spheroids, 
the thin one to the model by Chary et al. (2004, their Fig. 4), 
which includes the contributions from low luminosity starbursts 
(LiR < 10" L0, 3 dot-dashed), LIRGs with lO" < Lir < 10^^ 
L0 (dashed), and ULIRGs with Ljr > 10^^ Lq (dotted). 
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3.5 Contribution to Spitzer-IRAC counts 

In Fig. 1121 we show the contribution of spheroidal galaxies 
to the Spitzer/IRAC bands, as predicted by our reference 
model. As expected based on results for the K and the 6.7/im 
bands, spheroids dominate the counts (except perhaps at the 
bright end), with a contribution decreasing from the near- 
to the mid-IR. Spheroids in the active star-forming phase 
dominate at faint flux density levels ( ^ 30/iJy). 
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Figure 11. Expected distribution of the 8 /im flux density of 
our model spheroids brighter than 100 fi3y at 24 fim (reference 
model). 



3.4 Contribution to 24^m Spitzer counts 

Figure |5| shows the contribution from forming spheroids to 
the 24 jj,m number counts, and the predicted redshift distri- 
butions at various flux limits. 

In particular, our reference model predicts a non- 
negligible contribution, ~ 20%, to the 24 /xm Spitzer counts 
above ~ 0.11 mJy. Yan et al. (2004a) report that 13% of the 
detected sources do not have a counterpart at both 8 fj.m 
and 0.7 /im, to flux limits of 20 and 0.17 /xJy respectively. In 
our view, this sub-population should be dominated by form- 
ing spheroids. Figure lTTl reports the distribution of S,/(8/im) 



for our model population with S24p 



100 fiJy, and shows 



4 DISCUSSION AND CONCLUSIONS 

Granato et al. (2001, 2004a) have shown that the mutual 
feedback between star-forming spheroidal galaxies and the 
active nuclei growing in their cores can be a key ingredi- 
ent towards overcoming one of the main challenges facing 
the hierarchical clustering scenario for galaxy formation, i.e. 
the fact that the densities of massive high redshift galaxies 
detected by SCUBA and by deep near-IR surveys are well 
above the predictions. 

The model by GDS04 interprets the SCUBA galaxies 
as proto-spheroids in the process of forming most of their 
stars. In fact, according to this model, essentially all mas- 
sive halos virializing at z Z 1.5 give rise to elliptical galaxies 
or to bulges of later type galaxies, that form at lower red- 
shifts. The star-formation in small halos is suppressed by 
supernova-driven winds, so that the early star-formation oc- 
curs primarily in the more massive halos, through gigantic 
starbursts with star formation rates ~ 10^ Mq /yr, whose 
duration is determined by the feedback from the active nu- 
clei, followed by essentially passive evolution. The model 
can therefore be appropriately tested primarily against sur- 
veys selecting high- 2; galaxies, like the (sub)-mm surveys by 
SCUBA and MAMBO, or spheroidal galaxies, like near to 
mid-IR surveys. 

To this end we first assessed the values of the ad- 
justable parameter tmc of the spectrophotometric code 
GRASIL (Silva et al. 1998), which regulates the effect of 
the complex and poorly understood radiative transfer pro- 
cesses on the time-dependent SEDs. To properly deal with 
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the extreme conditions prevailing during the intense star- 
formation phase of spheroidal galaxies, a new parameter, k, 
had to be introduced. Having fixed the values of the two 
parameters, and keeping the GDS04 choice for the parame- 
ters controlling the star-formation history, the chemical en- 
richment and the evolution of the dust and gas content of 
massive spheroidal galaxies, we have shown that the model 
accounts, on one side, for the SCUBA and MAMBO counts, 
corresponding to the active star-forming phase, and, on the 
other side, for the if -band counts for K 20-21 and for the 
redshift distribution oi K < 20 galaxies, as well as for the 
high-z tails of redshift distributions up to if = 24. 

Although observations unambiguously indicate that 
stellar populations in nearby spheroidal galaxies are old, it 
is possible that stars were formed in small sub-units, be- 
fore the galaxy was assembled. Our results highlight a clear 
continuity between objects where stars formed (detected by 
(sub)-mm surveys) and the evolved galaxies dominating the 
bright if-band counts, indicating that massive spheroidal 
galaxies formed most of their stars when they were already 
assembled as single objects. As illustrated by Fig. |21 the 
model is remarkably successful in reproducing the observed 
K < 20 redshift distribution for z > 1, in contrast with 



both the classical monolithic models (which overestimate the 
density at high-z) and the semi-analytic models (that are 
systematically low). For example, Somerville et al. (2004) 
find that their hierarchical model underproduces K 20.15 
{Kab ^ 22 in their notation) galaxies by about a factor of 
3 at 2 ;j 1.7, and by an order of magnitude at z ^ 2, while 
a monolithic model overproduces these galaxies by about 
a factor of 2 at z ~ 2. The GDS04 model also reproduces 
the bimodal distribution of colours of the Somerville et al. 
(2004) sample, as due to contributions of star-forming and 
passively evolving spheroids. The ratio of star-forming to 
passively evolving spheroids is also nicely consistent with the 
finding of Somerville et al. (2004) and Daddi et al. (2004), al- 
though the observations of the latter authors suggest a more 
gradual decline of the star-formation rate than implied by 
the current version of our model. Again, reproducing the ob- 
served colour distribution proves to be very challenging for 
both monolithic and standard semi-analytic models. 

A specific analysis of the K-hand counts and redshift 
distributions of EROs has been worked out, since these 
objects appear to provide particularly sharp constraints 
on galaxy evolution models. The results confirm that the 
GDS04 model successfully reproduces the counts, and is in 
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reasonably good agreement with the redshift distributions, 
although the need of a more sophisticated modelling of the 
star formation history and of dust geometry is indicated. 

The mid-IR counts provided by ISO CAM, primarily 
at 6.7 and 15 /im, and by the Spitzer Space Telescope at 
several additional wavelengths and over larger areas, con- 
strain the SEDs of spheroidal galaxies which dominate for 
A ^ 6 fj,m, and the evolution properties (as well as the SEDs) 
of late-type (normal and starburst) galaxies that dominate 
at longer wavelengths. The ISOCAM counts at 6.7/im and 
the associated redshift distribution are well fitted by the 
model (Fig. |KJ, as are the Spitzer/IRAC counts at 3.6, 4.5, 
and 5.8 ^m. 

At longer wavelengths, we expect that the contribution 
of starburst galaxies to the counts becomes increasingly im- 
portant, particularly at bright flux densities, and to actually 
dominate at 24 /im (Fig.l^J- Although modelling the counts 
of such sources is beyond the scope of the present paper, 
we note that, according to the GDS04 model, essentially all 
massive halos at z Z 1.5 turn into spheroidal galaxies. As 
a consequence, we expect a sharp drop in the redshift dis- 
tribution of starburst galaxies at z <; 1.5, in contrast with 
predictions of some current phenomenological models (Dole 
et al. 2003; Chary & Elbaz 2001), and a high-z tail due 
to star-forming spheroids. For example, we expect a surface 
density of sources with 5*24^111 > 20 fiJy and z ~ 2 about 
2 times lower than predicted by Chary et al. (2004, their 
Fig. 4). 

While the scanty redshift information on 24 /im sources 
(Le Floc'h et al. 2004) does not allow us to test our expec- 
tations yet, such sources can be preliminarily characterized 
by their IR colours (Yan et al. 2004a). Remarkably, the pre- 
dicted fraction of star- forming spheroids with 8 fim fluxes 
below the detection limit matches that found by Yan et al. 
(2004a). We thus interpret these sources as proto-spheroidal 
galaxies in the process of forming most of their stars, at 
typical redshifts ~ 2 but with a tail extending up to z ~ 4. 

Moving to still longer wavelengths, we expect that star- 
forming spheroids start dominating the 160 /im counts below 
50-100 mjy, and produce a peak at z ^ 1.5 in the redshift 
distribution, while at brighter fluxes, counts are likely domi- 
nated by low redshifts sources (Rowan- Robinson et al. 2004) . 

It is worth noticing that the AGN activity powered by 
the growing SMBHs has a small effect on the number counts 
of spheroids at all wavelengths considered here. Since the 
AGN SED {uLv) is relatively flat at optical to far-IR wave- 
lengths and drops in the sub-mm region, the maximum AGN 
contribution occurs around 15 ^m, at the minimum between 
the peak of direct starlight emission and that of dust emis- 
sion (Silva, Maiolino & Granato 2004; see also Fig.|SJ. But, 
at this wavelength, the spheroids are anyway a minor compo- 
nent of the counts (see the upper left-hand panel of Fig. |3Ji. 
The sub-mm counts are dominated by proto-spheroids with 
very high SFRs (SFR ^ 500 Mq yr for sources brighter than 

5 mJy), while the accretion rate onto the SMBH predicted 
by the model is several orders of magnitude lower (see Figs. 2 
and 3 of GDS04). Therefore the AGN luminosity is generally 
much lower than that due to stars (it is non-negligible only 
for a very small fraction of the star forming phase) even al- 
lowing for a much higher mass to light conversion efficiency. 
On the other hand, we expect the mild AGN activity to be 



detectable in X-rays in the majority of the bright SCUBA 
sources (Granato et al. 2004b), in keeping with the find- 
ings by Alexander et al. (2003). A detailed analysis will be 
presented in a forthcoming paper (Granato et al. in prepa- 
ration) . The counts of spheroids at A < 10 yLxn are domi- 
nated by objects in the passive evolution phase, for which 
the mass accretion rates onto the SMBHs are very low and, 
correspondingly, the AGN activity is very weak. 
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